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Abstract 

Using Painleve singularity structure analysis, we show that coupled higher- 
order nonlinear Schrodinger (CHNLS) equations admit Painleve property. Using 
the results of Painleve analysis, we succeed in Hirota bilinearizing the CHNLS 
equations for the integrable cases. Solving the Hirota bilinear equations, one 
soliton and two soliton solutions are explicitly obtained. Lax pairs are explicitly 
constructed. 



PACS number (s): 42.81. Dp, 02.30 Jr, 04.30 Nk 



1) Electronic address: vinojmnOcusat.ac.in 

2) Electronic address: vck@cusat.ac.in 



(to appear in Phys. Rev.E) 



2 



I. INTRODUCTION 



The last three decades witnessed extensive theoretical and experimental studies on 
optical solitons because of their potential applications in long distance communica- 
tion. The invention of high-intensity lasers helped Mollenauer et al [1] to verify ex- 
perimentally the pioneering theoretical work on optical solitons initiated by Hasegawa 
and Tappert [2] . The solitons, localised- in-time optical pulses, evolve from a nonlinear 
change in the refractive index of the material, known as Kerr effect, induced by the light 
intensity-distribution. When the combined effect of the intensity-dependent refractive 
index nonlinearity and the frequency-dependent pulse dispersion exactly compensate 
each other, the pulse propagates without any change in its shape, being self-trapped by 
the waveguide nonlinearity. The propagation of optical solitons in a nonlinear disper- 
sive optical fibre is governed by the well-known nonlinear Schrodinger (NLS) equation 
of the form: 



where u is the complex amplitude of the pulse envelope, a and (3 are the group velocity 
dispersion (GVD) and self-phase modulation parameters respectively and subscripts z 
and t represent the spatial and temporal coordinates respectively. 

When ultrashort pulses (USPs) are transmitted through fibres, higher-order effects 
such as third-order dispersion (TOD), Kerr dispersion and stimulated Raman scattering 
(SRS) come into play as experimentally reported by Mitschke and Mollenauer [3]. The 
Kerr dispersion, also known as self-steepening, is caused by the intensity dependence 
of the group velocity which results in asymmetrical spectral broadening of the pulse 
since the peak of the pulse travels slower than the wings. The SRS causes a self- 
frequency shift which is a self-induced red shift in the pulse spectrum as the low- 
frequency components of the pulse obtain Raman gain at the expense of the high- 
frequency components. With the inclusion of all these effects, Kodama and Hasegawa 
[4] have proposed that the dynamics of femtosecond pulse propagation be governed 
by a higher-order NLS (HNLS) equation. The HNLS equation allows sohton-type 
propagation only for certain choices of parameters [5]. 

A coupled NLS equation was proposed by Manakov, by taking into account the 
fact that the total field comprises two fields with left and right polarizations [7]. The 
coupled equation takes the form: 



iut-\- a Uzz -\- li\u\ u = Q 
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The above equations are integrable only for the following parametric choices: (i) 
ci — C2-, a — (3 — ^ and (ii) ci = — C2, a — (3 — ^. Recently, for USPs, Eq.(2) was 

generalized to a set of coupled higher-order NLS (CHNLS)equation which can be de- 
rived from the Maxwell's equations in order to investigate the effects of birefringence on 
pulse propagation in femtosecond regime [8, 9]. The general form of CHNLS equations 
are: 
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In general, the above equations are not completely integrable. However, if some 
restrictions are imposed on the parametric values, one can obtain several integrable, 
soliton- possessing NLS-type equations: (i) A = 0, NLS; (ii) /3i : /32 : /Ss = : 1 : 1, 
derivative NLS [10]; (iii) /3i : /32 : /^a = : 1 : 0, derivative mixed NLS [10]; (iv) A : /32 : 
(33 — 1 : 6 : 0, the Hirota equation [11] and (v) /?i : /?2 : /?3 = 1 : 6 : 3, Sasa-Satsuma 
equation [8]. Along the lines of Refs.[12, 13] we choose Pi — 1, P2 — 6, P3 — 3. For this 
choice of parameters, Eqs. (3) become 
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The plan of this paper is as follows. In Sec. II, we establish the Painleve property 
of the above system of equations. In Sec. Ill, we rewrite Eqs. (7) in a Hirota bihnear 
form. Section IV is devoted to the construction of exact one soliton and two sohton 
solutions. In Section V, Lax pair for CHNLS equations are obtained. Conclusion and 
discussion of the present calculation are presented in Section VI. 



II. PAINLEVE ANALYSIS OF CHNLS EQUATIONS 

In this section we study Painleve analysis of Eqs. (4). The motivation behind this 
exercise is that P condition is necessary one for studying the integrablity of nonlinear 
partial differential equations [14-19] and helps construct solutions. The method for 
applying Painleve test to partial differential equations were introduced by Weiss, Tabor, 
and Carnavale [15] with simplification due to Kruskal [16] involves seeking a solution 
of a given partial differential equation in the form: 
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u{z, t) = r T,T=o Uj{t)(l>i{z, t),Uo^O 

(5) 

v{z,t) = ^^i:T=ovAtW{z,t),vo ^ 

with 

= z + = (6) 

where il^it) is an arbitrary analytic function of t, Uj{t) and Vj{t) ( j = 0,1,2,....,), 
in the neighbourhood of a noncharacteristic movable singularity manifold defined by 
= 0. 

Apart from providing the integrability property of a given nonlinear partial dif- 
ferential equations, P analysis also provides information about Backlund transforma- 
tion(BT), Lax pair, Hirota's bilinear representation, special and rational solutions, 
etc. [15-17]. Many of these results are obtained by truncating the Laurent series at a 
constant level term [18, 19]. 

In order to investigate the integrability properties of Eqs.(7), we rewrite it in terms 
of four complex functions a, b, c and d by defining u — a, u* — b, v — c, v* — d. 
Consequently, we have the following equations: 

iat + ttzz + '2{ah + cd)a — i\[azzz + 6(a6 -|- cd)az + ^{ab + cd)za\ = 
-ibt + bzz + 2(a6 + cd)b + iX[bzzz + 6(a6 + cd)bz + 3(a6 + cd)zb] = 

(7) 

ict + Czz + 2{ab -\- cd)c — iX[czzz + 6(a6 -|- cd)cz + 3{ab + cd)zc] — 

—idt + dzz + 2(a6 + cd)d + iX[dzzz + 6(a6 -|- cd)dz + 3(a6 + cd)zd] = 

The Painlcve analysis essentially consists of four main stages (i) determination of 
leading-order behaviour, (ii) identifying the resonance values, (iii) verifying that at 
resonance values sufficient number of arbitrary functions exist without the introduc- 
tion of movable critical manifold, and (iv) identifies connection with the integrability 
properties hke Lax pair, BT. 

Looking at the leading order behaviour, we substitute a ~ ao0°\ b ^ bo0"^) c ~ 
Co0"^, d ~ (^00"" in Eqs.(7) and balancing the different terms, we obtain the following 
results: 

(Xi — a2 — 013 — a4 — —1, 

(8) 

aobo + codo = 
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For finding the powers at which the arbitrary functions can enter into the series, we 
substitute the expressions, 

a — aQ(f)~^ + aj(p^^^, b = bo(p~^ + bj(/)^~^ 

(9) 

c — Co(j) ^ -\- Cj(j^ ^jd — docf) ^ + dj(f)^ ^ 

into Eqs. (7), and compairing the lowest-order terms we obtain a system of four hnear 
algebraic equations in {aj , bj , cj , dj ). In matrix form it may be conveniently written 
as 

[A{j)][X] = (10) 



where [X] 
and 



[Am = 



[ttj , bj , Cj , dj ) 
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C 

-6dl 
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where 

B^U- l)(j - 2)(i - 3) - 3{j - 1) - 6ao6o + 3(i - 2) + 3 and 
C=ij- l)(j - 2)(i - 3) - 3{j - 1) - Gcodo + 3{j - 2) + 3 
To have a non-trivial solution for aj , bj , Cj and dj,we demand that 



+3aoCo{j - 2) 

-66oCo 
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C 



det A{j) = 



On solving Eqs. (11), we get the resonance values ets j — —1, 0, 0, 0, 2, 2, 2, 3, 4, 4, 4, 4. 
The resonance at j = -1 corresponds to the arbitrariness of ip {z, t) . On equating the 
coefficients of i^""^, we get a unique equation defining oq, &o,Co and do which is given by 

aobo + codo ^ ~. (12) 

This shows that any three of the four functions ao, &o,co and do are arbitrary which 

corresponds to j = 0, 0, 0. 

Proceeding further and equating the coefficients of 4>~^ , V'"^) , we obtain 
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"1 3iX ' 
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On the other hand, the coefficients of {ip ^,'0 '^,ip "^^il) ^) in Eqs. (7) reduce to a 
single equation 



oA 



(14) 



so that three of the four functions 02, 62, C2 and (i2 are arbitrary which corresponds to 
J = 2,2,2. _ _ _ _ 

Similarly from the powers of {-^'^ ^-^'^ ,il)~^)aindi (-0"°, V'"^) 1 we find 
that Eqs. (4) admit the sufficient number of arbitrary functions and hence Eqs. (7) 
possess the Painleve property and hence they are expected to be integrable. 

III. HIROTA BILINEARIZATION 

Hirota's bilinear method [20] is one of the most direct and elegant methods available 
to generate multi-soliton solutions of nonlinear partial differential equations. To avoid 
mathematical complexities, it is rather convenient to transform Eqs. [4] to a simpler 
form, so that we may able obtain multi-soliton solutions. We make the following 
transformations to convert CHNLS to complex modified K-dV (cmK-dV) equation: 



u {z, t) = Qi {Z, T) Exp [-1 (f - ^)\ , V {z, t) = Q2 (Z, T) Exp [-^ (f - ^) 



t = T,Z = z + l, 



(15) 



Using the above transformations in Eqs. [4], the resultant cmK-dV equation is obtained 
in the form 



QiT - A [Qizzz + 6 (iQif + IQ2I') Qiz + 3gi {\Qif + \Q2f) 
Q2T - A [Q2ZZZ + 6 (iQil' + IQ2I') Q2Z + {\Qif + \Q2 



= 







(16) 



In order to construct Hirota's bilinear form, we consider Hirota bilinear transformations 
in the form 
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Qi = - Q2 = (17) 

where G {Z,T) and H {Z,T) are complex functions and F {Z,T) is a real function. 
Now using the transformations (17), (16) can be rewritten as 

[(Dt - \D%) {G.F)\ - \[{-?>Dl (F.F) Dz + 12 (\Gf + \Hf) Dz) (G.F) + 
SGFDz (G.G*) + 3H*FDz {H.G) - SHFDz {G.H*)] = 

(18) 

[{Dt - XD%) (H.F)] - X[{-3Dl (F.F) Dz + 12 (|G|' + \Hf) Dz} (H.F) + 
SHFDz (H.H*) - 3G*FDz (H.G) + SGFDz {H.G*) = 
where the Hirota bilinear operators D^ and Dt are defined as 



D^D-G {Z, T) F {Z', T') = [q^ ' g^j [of ' Qr) ^ {Z, T) F {Z' , T') \z=z',t=t' 

(19) 
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IV. EXACT SOLITON SOLUTIONS 

Eqs.(18) can be decoupled into a set of bilinear equations as: 

(Dt - \D%) {G.F) =0 {Dt- \D%) {H.F) = 
D|(F.F) = 4(|G|' + 
Dz (G.G*) = 0, Dz (H.H*) = 0, Dz (G.H*) = 0, Dz (H.G*) = 0, Dz (G.H) = 

(20) 

In order to obtain soliton solutions, we arc applying a perturbative technique by writing 
the variables F, G, H as a. series in an arbitrary parameter e: 

F = 1+£V2+£V4+ , G = egi+e^g3+e^g5+ , H - ehi+e^h3+e^h5+ 

(21) 

A. One-soliton solutions 

For one-soliton solution (ISS), we assume solutions in a series form in s such that 

F = 1 + 6^2 G^egi H ^ ehi (22) 

Substituting Eq.(22) in Eqs.(20) and then collecting coefRcents of terms with same 
powers in e, we obtain: 
e : 

(Dt - XD%) (^1.1) = (Dt- \D%) (/ii.l) = (23) 



£2 



Dl{l.f2 + f2.l)^^{gi.gl + hi.hl) 

(24) 

Dzigi.gl) = 0,Dz{h.hl) = 0,Dz{gi.hl) = Q,Dz{h^.gl) = 0,Dz{gi.h) = 



[Dt - \Dl) {g^.f^) = {Dt- \Dl) {h^.h) (25) 



Dl (/2-/2) (26) 
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One can easily check that the solution, which is consistent with the system (23-26), is 

gi = cos exp (77) hi = sin exp (77) 

h - exp {2rj) 

where 



(27) 



r] = kZ + Xk^T 



(28) 



and (f) and k are real constants. Using Eqs. (27) and (28) in (22) and then in (17), 
after absorbing e the one-soliton solution can easily be worked out to be 



= cos 0sec/i {kZ + Xk^T + rjo) 
Q2 = (^) sin 0sec/t {kZ + Xk^T + r]o) 



(29) 



where r]o is a real constant. Using Eqs. (15), one-soliton solutions of Eqs. (4) are found 
to be 



u 



{,72} "^^^ ^ ^-^P ~^ ~ 2^) (^^ + -^^^^ + ^0) 
sin exp [-i (^|^ - 2^)] sec/i (/cZ + Xk^T + rjo) 



(30) 
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B. Two-soliton solutions 



The two-soliton solutions can be obtained by terminating the series expansion for 
F, G, H as: 



(31) 



and proceeding as before to obtain 

gi = exp (r/i) + exp (7/2) h = igi 



93 = {h - hy 



exp{2rii +ri2) i exp(r]i +2??2 ) 



f — A \T^p(23i1 _i_ o exp(??i+??2) I exp(2??2) 



(32) 



where 



/4 = 



4(fc2-fci)'^exp(2?7i+2?72) 
16fe2fe2(jfcj+fe2)4 



r]j ^ kjZ + Xk^jT J = 1,2 



(33) 



Here A;^ is a real constant. Using Eqs.(32) and (33) in (31) and then in (17), the 
two-soliton solutions of (16) are obtained. Using Eqs.(15), the two-soliton solutions of 
Eqs.(4) are found to be: 



u 



f exp 



V — y exp 



Hi 



27A2 ; 



(34) 



Both ISS and 2SS are in exact agreement with Eqs.(4).Two dimensional and three 
dimensional plots of 2SS are given in figs.(l) and (2) respectively. 
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V. LAX PAIR FOR CHNLS SYSTEM 

The hnear eigen value problem associated with Eqs.(16) are [21,22] 



(35) 



where 



U ^ 



( -ix Qi Q*i Q2 Q*2 \ 
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-Qi 

V -Q2 
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/ ai2 ai3 ai4 015 \ 
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(37) 



where 

ai2 = 4 A Qi + Q^zz, ai3 = 4 A Q\ + Q\zz. = 4: A Q2 + Q2ZZ, 
ai5 = 4 A + Qlzz 

(I2I — —4 A QI — Q\zz^ ^22 = QlQlZ ~ QlQlZ, 0^24 — Q2Q1Z ~ Q1Q2Z, 0,25 = Q2Q1Z 
Q1Q2Z 

^31 = ~4 A Qi — Qizz, 033 — QiQiz — QiQizi ^34 = Q2Q1Z — Q1Q2Z, 035 — Q2Qiz 
Q1Q2Z 
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<^41 — ~^ ^ Q2 ~ Q2ZZ1 ^42 = Q1Q2Z ~ Q2Q1Z, Q43 — Q1Q2Z ~ Q\Q\zi ^^44 = Q2Q2Z ~ 

^51 = —4 AQ2 — Q2ZZ1 O52 = Q1Q2Z — Q2Q1Z, 0,53 — Q1Q2Z — <^55 = Q2Q2Z — 

Q2Q2Z 

with A = \Qif + \Q2\\ 

The compatibility condition Ut — Vz + [U,V] — gives rise to Eqs.(16). The 
construction of Lax pair confirms that Eqs.(16) and thereby the CHNLS Eqs.(4) are 
indeed completely integrable. 

VI.CONCLUSION 

In this paper, we have considered a set of coupled NLS equations with higher-order 
linear and nonlinear dispersion terms included. Then, by choosing the parameters as 
in the case of the corresponding integrable uncoupled case, we applied the Painleve 
singularity structure analysis and established that for this particular choice of param- 
eters, Eqs.(4) possess the Painleve property. We have explictly obtained one-sohton 
and two-soliton solutions for the integrable cases of CHNLS equations using Hirota bi- 
linearisation technique and solutions are plotted. We have also constructed Lax pairs 
using AKNS formalism. Hence, with these results, we have proved that the CHNLS 
equations which describe the wave propagation of two fields in fibre systems with all 
higher-order effects such as TOD, Kerr dispersion and stimulated Raman effect, will 
allow soliton-type propagation. From the soliton solutions, one can obtain information 
about the shape, width and intensity of the propagation pulse. 
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Figure captions. 

1. Figure la : 3D profile of \u{z, t) \ for the two-soliton solution of Eq. (34) with the 
parameter values ki=0.034, k2=0.04, A = 0.005 

2. Figure lb : 2D profile of \u{z, t) \ for the two-soliton solution of Eq. (34) with the 
parameter values ki=2, k2=3, A — 0.002 

3. Figure Ic : Contour plot of |it(z, t)| with respect to z and t for the parameter 
values ki=0.034, k2=0.04, A = 0.005 

4. Figure 2a : 3D profile of \v{z, t) \ for the two-soliton solution of Eq. (34) with the 
parameter values ki=0.04, k2=0.045, A = 0.005 

5. Figure 2b : 2D profile of \v{z, t) \ for the two-soliton solution of Eq. (34) with the 
parameter values ki=2, k2=3, A = 0.005 

6. Figure 2c : Contour plot of |M(z,t)| with respect to z and t for the parameter 
values ki=0.04, k2=0.045, A = 0.005 
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